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As often stated, a class of f(R) gravity models seem to obtain ACDM behaviour for high red-shifts, 
and general relativistic behaviour locally at high curvatures. In this paper, we study numerically 
sun-like configurations for one such proposed model, namely Hu-Sawicki f(R) model having the 
chameleon property. 

Chameleon gravity theories are stated to reproduce general relativistic behavior within high cur- 
vature regime. Therefore stationary spherically symmetric configurations of f(R) gravity with poly- 
tropic equation of state and regular metrics at the origin have been studied numerically and com- 
pared with general relativistic Eddington polytropes. 

To determine the metric deviation from relativistic Schwarzschild-de Sitter solution, post- 
Newtonian parametrization is used around the f(R) polytropes. Since, Birkhoff 's theorem does not 
need to hold within a modified gravity theory, the solution outside may be different. We find that 
although Schwarzchildian space-times can be produced with this model, they demand fine-tuning. 



I. INTRODUCTION 

According to observations [TJ [2] , a dark energy (DE) 
component is needed to explain the apparent acceleration 
of the Universe. In the ACDM model, Einstein's general 
relativity (GR) is modified by a tiny constant A. With 
this deviation from the bare Ricci scalar in the Einstein- 
Hilbert action, the universe accelerates according to cos- 
mic background, supernovae and large scale observations 
[DS]. Various other means to introduce the accelerated 
expansion into the field equations include e.g. cosmo- 
logical scalar fields with gravitational interactions [5], as 
well as vector fields [6], string theory and higher dimen- 
sion inspired models [7] as well as several type of gravity 
modifigations [8]. 

f(R) modified gravity theories [HI UH], where the 
Einstein-Hilbert action is generalized as a function of the 
Ricci scalar, provide natural alternatives for dark energy 
that needs new physics to explain the accelerated expan- 
sion of the universe. In these models some subdominant 
terms in the action functional become essential in very 
large scales. 

In the present paper, the f(R) theories of modified 
gravity are studied, that include an action function with 
a suitable function of the Ricci scalar. The observed 
large scale defects could then be explained by the devi- 
ation from Einstein-Hilbert gravity. However, Einstein- 
Hilbert gravity is well tested within solar system length 
scales and at low energies [TT]. Therefore we, focus on 
models considered as alternatives for the concordance 
ACDM model with GR like behaviour at solar system 
scale. 

When the modification in the action functional is writ- 
ten using auxiliary scalar fields (as in quintessence and 
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scalar-tensor teories [8]), the new effects of gravity 
can be interpreted as scalar-field - matter interactions. 
Chameleon f(R) models [T5] allow the mass of the scalar 
field to vary according to the local matter density and 
give rise to non-linear self-interactions. With chameleon 
behaviour, correction to the Newtonian gravitational po- 
tential may be small in high curvature regimes (e.g. near 
the sun) and the metric therefore close to Schwarzschild- 
de Sitter solution. 

There are many f(R) models that can produce correct 
shift from matter domination to radiation domination 
and to the currently accelerating epoch (see e.g. p~3 | I14 J ) . 
Also, constraints from structure formation and stability 
conditions are satisfied by various models [T3HIS]- Since 
gravity is most rigorously tested at local scales and with 
low energies, tangible f(R) model must be able to set- 
tle to the observed GR values within this regime. f(R) 
models that pass all the above mentioned tests are scarce. 
Among current viable f(R) models are [20Tf24] . In the 
present paper we study Hu and Sawicki /(i?)-model (HS) 
numerically to test if it locally deviates from GR. It has 
non-linear rational action functional and seems to follow 
locally the chameleon mechanism. 

While the metric outside a spherically symmetric body 
in GR must be the Schwarzchild-de Sitter space-time ac- 
cording to Birkhoff 's theorem [35], this does not hold 
for a spherically symmetric matter distributions in non- 
trivial f(R) theories, like in the HS-model. In these the- 
ories, the way space-time outside the mass distribution 
depends on the solution inside, is more complex [57] . 
Therefore, we have compared static spherically symmet- 
ric bodies with a polytropic equation of state within the 
HS model and in Einsteinian gravity. 

We also calculate the mass of auxiliary scalar degree of 
freedom and find it to be high enough to evade the "fifth 
force" problem [T3] , 

This paper is organized as follows: first in Section [TT 
the gravitational framework is discussed, in Section [TTT 
the polytropic configuration and the surrounding space- 
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time is constrained and in Section Hvl the numerical work 
is represented. Finally in Section [V] we discuss the re- 
sults. 



II. GRAVITY FORMALISM 

We consider the f(R) gravity action [20] of the form 

s = J d A x + c m ) , (i) 

where k 2 = 8irG, and C m is the matter Lagrangian. The 
corresponding field equations, derived in the metric for- 
malism, are 

F(R)R M „ - ^f(R)g^ - (V^V„ - g M „n)F(R) = SirGT^. 

(2) 

Here is the standard minimally coupled stress-energy 
tensor for perfect fluid and F(R) = df/dR. Contracting 
the field equations we get the trace equation 

F(R)R ~ 2/(7?) + 30F(R) = 8ttG{ P - 3p). (3) 

which can be used as one independent dynamical equa- 
tion. 

As usual, the matter term needs to obey the equation 
of continuity D^T^" = 0. The only non-trivial compo- 
nent for spherically symmetric system is 



B' 



(4) 



where comma stands for radial derivative, ' = d/dr. 

Like in GR, the equation of continuity is automati- 
cally satisfied here if the field equations are satisfied [28] 
and no additional information is gained. In this work 
we chose the continuity equation to be one of the solved 
dynamical equations instead of using the full set of field 
equations. Therefore, as the set of independent equa- 
tions we use the 11-component of equation ([2]), the trace 
equation ([3]), the definition for the Ricci scalar in terms 
of metric components and the continuity equation Q. 



A. The Hu-Sawicki f(R) model 

As there are still several physically acceptable f(R) 
models we chose the Hu-Sawicki model for numerical pur- 
poses. This model is able to reproduce ACDM at high 
curvature when the f(R) function is expanded. It also 
reproduces correct matter domination and the current 
accelerated eras and avoids Ostrogradski and matter in- 
stabilities [T5HT7] [20] . This particular model also van- 
ishes for flat space-time, /(0) = 0, and is able to produce 
the asymptotic limit f(R)~R + A&sR— >oo. 

The action function for the Hu-Sawicki model is given 
by [20] 



f(R) H s = R - mf 



ci( 



R \n 



This contains three essential parameters, a positive real 
number ci, C2 and a positive integer n. The mass scale 
m can be always included in parameters c%, c 2 , but it 
is convenient to use it as the mass scale as in the paper 
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n m h 2 

0.13 



(6) 



where h is the reduced Hubble parameter and h = 0.72 
in our work. 

The authors of [2D] study the solar system behaviour of 
the model by an approximated f(R) function. They con- 
clude, that the model should correctly arrive at the solar 
system value outside a sun like matter distribution be- 
cause the scalar degree of freedom is locked to the general 
relativistic value out through the solar corona {20] . The 
approximate solution for ^ in high curvature (i.e. when 
m 2 /R — > 0) gives accordingly an effective ACDM model 
when the f(R) function is expanded. By approximating 
the f(R) function to begin with, the non-linearity of the 
modified equations that allow the Birkhoff's theorem to 
be violated, is not taken fully into account. 

In the present paper the above assumption is tested 
by using the full non-linear, rational action function with 
the non-linear field equations. We, therefore, numerically 
study whether the approximation of the rational function 
([5 ) , that makes the solution fully merge with the general 
re ativistic space-time around a polytrope, is appropri- 
ate. This is important because the model is currently 
used as a viable DE model (see e.g. [2"D]1. 

One useful way to study this family of f(R) models is 
to write ^ as: 



f(R) = R 



C2 



cim 



C 2 1 + C2 



(7) 



(5) 



Now one finds that only when the Ricci scalar R is near 
the vacuum value m 2 , can the third term contribute to 
ACDM . However, because of the non-linearity of the field 
equations, behaviour of the solution may be different 
from ACDM . 

The system was first required to give the correct static 
vacuum value R vac = 12Hq with the trace equation 
This restricted the initial-parameter space considerably. 
Moreover only analytical and regular solutions at the ori- 
gin were considered, by demanding the solutions to reach 
a finite value at the origin [3D]. 



III. STATIC SPHERICALLY SYMMETRIC 
SOLUTIONS 



We consider spherically symmetric, static configura- 
tions and adopt the metric for this configuration: 

ds 2 = B(r)dt 2 ~ A(r)dr 2 - r 2 (d6 2 - sin 2 6 d(f> 2 ). (8) 

SdS-metric describes an expanding vacuum solution of 
the Einstein field equations. In order to test whether 
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the studied modified gravity polytrope can really repro- 
duce the general relativity weak field limit, the metric of 
the HS-polytrope is compared to the SdS-metric on the 
boundary of the spherical matter distribution, 

Further, parametrized Post-Newtonian parameters are 
compared to observations in order to determine the rele- 
vance of the SdS-solutions in the solar system. Derivation 
of the 7pp n parameter by fitting the general metric ([8]) 
to the PPN SdS-solution, was done in [3U]. The used 
parameter is 



1 + 



iPPN = 



B(r) 
rfl'(r) 



1 



Mr) 
rA'(r) 



(9) 



here the cosmological term H is neglected because it is 
vanishingly small within solar system. 



A. Polytropic considerations 

We use standard polytropic configurations that resem- 
ble young sun-like solar mass main sequence stars. We 
chose a solar polytrope as the reference star with poly- 
tropic index n p — 3. This is the so called Eddington 
model and it is given by the general relativistic Tolman- 
Oppenheimer-Volkoff (TOV) equation. This model is 
able to reproduce the observed solar mass and radius 
rather well with a simple set of equations and initial con- 
ditions. 

As self gravitating globes of gas, the stars are kept 
in hydrodynamic equilibrium with the equilibrium con- 
dition 



dp 
dx a 



= P 



F a . 



(10) 



The TOV equa- 



Here F a is the gravitational force [ 
tion is defined with the hydrostatic equilibrium ( 10 ) and 
spherically symmetric metric ([8]). 

As static stellar objects, we consider perfect fluid solu- 
tions with polytropic equation of state and parametriza- 
tion 



P 
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(11) 
(12) 



The index n p = 3 gives the Eddington polytropic solution 
for the Lane-Emden (LE) equations for solar mass M 
and radius Rq. 

According to [31] the solar polytropic model with in- 
dex n p = 3 is not stabile with respect to Jacobi stability 
analysis. The solar polytropic index is stabile within lin- 
ear stability analysis and Kosambi-Cartan-Chern theory 
(see [31] and references therein). The polytropic index 
values: 



(11 + 8\/2)/7 < n p < 5 



Mission 


Measured PPN-parameter constraints 


Cassini 


7 - 1 = (2.1 ± 2.3) x 10" 5 


VLBI 


7-1 = (-2 ±4) x 1CT 4 


LLR 


4/3 - 7 - 3 = (4.4 ± 4.5) x 1(T 4 


Mercury perihelion 


|2 7 -/3-l| < 3x 1(T 3 



TABLE I: Observational constraints from solar system mis- 
sions: Cassini spacecraft mission [33], Very Long Baseline In- 
terferometry (VLBI) 34 , Lunar Laser Ranging (LLR) tests of 
general relativity [35] and the Mercury perihelion shift mea- 
surements [361 [37]. The constraints for Pppn are not dis- 
cussed in this paper, as the Cassini mission results give the 
most stringent limits for SdS here even for /3ppjv = 1 



would conform to also Jacobi stabile polytropes. How- 
ever, when departing from the solar index the radius- 
mass relation changes dramatically. The behaviour of 
the polytropic star with the used field equations would 
also change for this index and new degrees of freedom 
would be needed in the analysis. 

The interconnected Lane-Emden coefficients (K, 7 in 

1 need to be 



(111 and scale length I in LE equations [! 
chosen correctly in order to obtain a star that best de- 
scribes the sun. Lane-Emden equations for a polytropic 
configuration are derived from the general relativistic 
Tolmann-Oppenheimer-Volkoff equations. We also de- 
rived general relativistic ACDM model, for reference, us- 
ing the ansaz described in Section ([II]). This model is of 
course identical to the LE polytrope. The standardized 
Eddington model was chosen for this study, because it is 
considered to be in good resemblance with the standard 
solar model [32] with the central density 



Pi 



Pq = 76.5 g cm 



(14) 



We varied the central density of the Eddington polytrope 
in order to find statistically significant differences or simi- 
larities between the modified gravity polytropes and gen- 
eral relativistic ones. Also, other realistic main sequence 
stars (Mhs 7^ Mq) were searched from larger polytropic 
parameter space with f(R) gravity. But as in the case of 
general relativistic n p — 3 polytropes, the only possible 
configurations are those with the solar mass. This occurs 
because (see e.g. (11.3.18) in [25] ) 



M 



(15) 



(13) 



fl(37-4)/( 7 -2)_ 

The mass M is constant for the Eddington model n p 
B. Space-time constraints 



If the studied f(R) modification (5j) is to produce gen- 
eral relativistic gravitation within the weak field of solar 
system, the Schwarzschild-de Sitter space-time is antici- 
pated as the resulting space-time around the calculated 
matter distributions. This metric is also a fair approx- 
imation around slowly rotating bodies, like the sun, so 
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we calculate the Post-Newtonian jp pn parameter at the 
boundary to see how well the general relativistic value, 
Ippn = 1, is reproduced. 

The Post-Newtonian parametrization (PPN) gives the 
strongest bound for a modified gravity theory within the 
solar system. Observational constraints are presented in 
the Table |!J The Cassini mission (331 137] measurement 
constrains the space-time parameter 7ppArmost strin- 
gently within the solar system: 



observational condition for jppn ( 16 ) when compared to 



Ippn - 1 = (2.1 ±2.3) x 10 



-5 



(16) 



Therefore, this Cassini constraint is used as reference for 
choosing the acceptable solutions that mimic general rel- 
ativistic ACDMpolytropes. 



IV. NUMERICAL STUDY 

We study the local effects of the f(R) model by inte- 
grating the matter distribution, p(r) from the center of 
the polytrope to the boundary. In the numerical integra- 
tion the center is defined to be under one meter. SI units 
with c = 1 were used throughout this study. The regular- 
ity at the origin was separately tested for the considered 
set of equations. 

The calculations were highly dependent on the equa- 
tion ansatz, thus we used one of the modified field equa- 
tions ([2| , equation for the Ricci scalar and the continuity 
equation Q to overcome numerical difficulties. 

We solved the dynamical equations for the metric 
components A(r),B(r), the energy density p(r) and the 
scalar curvature R(r), with respect to the radial coordi- 
nate r (as was done in [2l)l l3T)] ). 

The space-time around the polytrope is found by in- 
tegrating from non-singular central region with the most 
general spherically symmetric metric ^ outwards until 
a boundary is found at p(r) = 10~ 5 kg m~ 3 , as sun's 
photosphere lies at p ps = 2 x 10~ 4 kg m~ 3 . 

This definition gives a discontinuous surface, since we 
assume the space-time to settle to the vacuum SdS value 
at the defined surface of a polytropic star. The Solar 
system potential or galactic potential is not discussed in 
this paper, because all the observations here are made 
either inside the polytrope or at the surface. For ref- 
erence, the same set of equations in [TT] were also used 
with f(R) = R + A with the same numerical parameters. 
This is the Eddington polytrope and is also denoted as 
ACDM model. 



We parameterized the polytropes (111 as is done when 
deriving the Lane-Emden equations: 



Po 



flV(7- 
7 



(17) 



with (14 1, and could therefore variate around the solar 
polytrope parameters, i.e. pi and 7 were determined by 
the Eddington model that fits the standard solar model 
well and 9 was varied. This was done in order to find 
out how well the f(R) polytropes satisfy the stringent 



the Eddington polytrope. By variating #0 we were able 
to easily find statistically significant differences between 
these models, because "fpp n parameter depends on the 
central density. 

The model depends effectively only on the central den- 
sity 6*o in (17) as is shown in Chapter IV A since the 



scalar curvature effectively follows the energy density 



R(r) r~j np(r). 



(18) 



Similar parameterization was also done with the 
ACDM model and the resulting polytropes, although 
very similar, differ systematically. 



A. Constraining the parameter space 

For studying the possible space-times around the 
f(R) polytropes we scanned a large parameter space of 
(n, Ci,C2,Ro,0q). The initial values A$, B for the spher- 
ically symmetric metric coefficients (|8| were chosen by 
using the regularity condition at the center when r — > 0. 
The metric parameters 9 and </> ffify do not enter the equa- 
tion ansatz. 

We first constrained the parameter space of the f(R) 
function by considering only those parameters c\ and C2 , 
that give the correct static vacuum value for the Ricci 
scalar R = 12Hq. The allowed values include only solu- 
tions that reproduce the correct vacuum with the trace 
equation ([3| for 



Tuu = « p\. 



(19) 



The Hu-Sawicki parameters c\ and , due to this require- 
ment, have a fixed ratio. The parameters n,c\,c% were 
scanned over values n G [1, 12], c\ and C2 € [1, 10 10 ] and 
the results seem highly degenerate over these parameter 
values. 

00 was varied, in order to gain statistical significance, 
through the values that produce physical central densities 
of red giants po/po,Edd — 0.01 to a very dense star with 
central density Po/Po,Edd — 1000. 

The Hu-Sawicki model has also been compared to the 
Type la supernova (Snela) and gamma ray burst (GRB) 
data, by [38]. To note, the parameter values favoured by 
Snela and GRB data with 95% CL in that work, do not 
conform to the correct static vacuum solutions for this 
model. 



B. Behaviour of HS polytropes 

For the scanned range of initial scalar curvature val- 
ues the space-time outside the sphere is near to, but 
not exactly, either SdS or Einsteinian (corresponding to 
f(R) = R) space-time This is easy to see if the function 
f(R) is written as in Q where the last term vanishes for 



all initial values except near the vacuum R : 
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S(R), 9„ = 1.54 
Lane-Emden solution 

f(R), 6 Q = 0.6 




Radius jR . 



FIG. 1: General relativistic LE polytrope is compared to 
Hu-Sawicki f(R) polytropes. Radius of the /(i?)-polytropes 
changes inversely proportional to the central density param- 
eter 8q in (17 1, such that the mass of the polytrope is always 
the solar mass Mq. 



The density parameter gives the solar 
value with Go = 1. There also exists f(R) jjs-polytropes al- 
most identical to GR LE solution with yppN within observa- 
tionally acceptable limits (161. The /(i?)-polytropes shown 



in this figure, however, have the jppn value out of this ob- 
servable range and are not observationally favorable. 



Hu-Sawicki polytropic solutions are, indeed, very simi- 
lar to general relativistic polytropes. The energy density 
always tracks to the general relativistic ACDM solution 
throughout most of the stellar radius and, therefore, the 
radius and mass of the HS-polytrope are similar to the 
GR solution. 

A key issue that distinguishes GR and f(R)HS poly- 
tropes, is the space-time outside the sphere. The 
Ippn parameter is not sufficiently close to 1 around the 
f(R) polytropes for the majority of the solutions. 

The solutions show unstabile numerical behaviour for 
the scalar curvature R in the central regions and near the 
surface. The solution, however tracks to fulfill R(r) ~ 
Kp(r) for almost all the radius. Because of the tracking 
behaviour, initial scalar curvature can be almost any- 
thing and still produce a GR-like solution. Generally 
with a non-SdS space-time, the only effective initial pa- 
rameter for distinguishing Hu-Sawicki model from the 
ACDM model is, therefore, the central density parame- 
For the central density of the Eddington 
= 76.5 g cm -3 ), the space-time parameter 



ter O in ([17]). 
polytrope Jp Q 

outside the polytrope varies around the SdS-value with 
a mean value close to the ACDM value (see FIGJ2]). 

Because the equations in[H]are of higher order than the 
field equations of general relativity, the conditions at the 
surface of the star need to be more stringent. For exam- 
ple, the metric parameters B and A in (|8j), need not be 
continuous at the surface of the star in general relativity 



but in the case of fourth order gravity, are required to be 
continuous. The reference value for a general relativistic 
Eddington polytrope at the defined surface approaching 
from inside is jppn.acdm — 0.9999991717. It is exactly 
one once the vacuum is reached, since the polytropes sur- 
face is discontinuous at the border by the definition of our 
polytropes. 



C. Numerical results 

For a wide range of Hu-Sawicki parameters and initial 
values for the energy density parameter 8 in (171, the 



solution is similar to the general relativistic polytropes. 
The notable difference is the space-time around the con- 
figurations, that turns out to be non-SdS for roughly half 
the f(R) stars per scan. 

All the polytropic modeling was done with the full 
non-linear field equations and non-linear rational func- 
tion ([5]). The considered numerical accuracy did not suf- 
fice for finding any solutions for the used set of equations 
( |2|3|4| ) with n > 9. 

Overall, for n = 1 ... 8 in ^ with the chosen ansatz, 
higher than 13 decimal-accuracy was not feasible for a 
large set of solutions. This accuracy, however is suffi- 
cient for finding an adequate 7ppjv value. With higher 
accuracy most of the equations were not solvable because 
of the high degree of non-linearity. When the accuracy 
is improved, the amount of successful polytropes drops 
dramatically because of numerical difficulties. Even with 
increasing accuracy the mean 7pp^ value does not con- 
verge to the ACDM value, as is shown in FIGjSJ 

The solutions also have numerical instabilities in the 
Ricci scalar, R(r), near the center and surface. The so- 
lutions, however, always track to fulfill R(r) = np(r) ef- 
fectively as the resulting configuration. 

Outside the correct parameter area confronting the 
static vacuum solution for the trace equation, f(R) poly- 
tropes do not yield 7ppjv = 1/2 outside, as in the case of 
the simpler polynomial f(R) models. Rather, the jppn 
value lies near 1 everywhere except if the initial value for 
the scalar curvature is that of the vacuum. 

All the solutions were required to be regular at the 
origin. The regularity condition was not, however, suffi- 
cient to force the solution to be the Schwarzchild-de Sitter 
space-time This is because around 60% of the polytropes 
do not fulfill the observational requirement of SdS given 
by pi. 



The Hu-Sawicki parameters are degenerate when con- 
sidering the mass, radius and jppn of the polytrope. 
This value is, however, dependent on the central density 
of the polytrope. To get statistically relevant comparison 
to the ACDM model with constant 7pp n , polytropic so- 
lutions with modified gravity were scanned over a wide 
range of central densities. 

As a result, the PPN parameter values typically lie 
within |7ppat — 1| < 1 x 10 4 , e.g. the solutions usu- 
ally tend to a SdS-like space-time. SdS space-time is, 
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FIG. 2: ("/ppn -1) value is plotted with respect to varied cen- 
tral density 8o in ( |17[ ). Only the observationally acceptable 
values according to Cassini mission (161, and around the so- 
lar polytrope, are shown. Lower curve shows the mean value 
<"fppN ~ 1 > per initial central density 9q. For comparison, 
origo line and a constant line (the upper white line) showing 
the ACDM model (7ppjv-1) value that was also calculated 
with our code. 



however, rarely reached. With only 10%, in a typical 
29000 polytrope run, converging to (jppn — 1) < 10~ 6 . 
With a mean value (jppn — 1) = — 6 x 10~ 6 . The 
ACDM value calculated with f(R) = R + A in the code 
gives (jppn^cdm — 1) = — 8 x 10~ 7 . Furthermore, only 
one percent of the HS solutions have (jppn — 1) < 10~ 7 . 
Typically 40% of the solar (i.e. 9q ~ 1) solutions at the 
surface (and 30% of all the solutions) have the SdS space- 
time parameter within the acceptable observational lim- 



its (161 



In FIG|2] mean values for the space-time parameter 
Ippn is calculated over varying initial central density 9 . 
Only the observationally acceptable values are plotted. 
In FIG. [3] the whole set of polytropes is shown. In both 
figures, the ACDM value is also plotted for reference. It 
is always closer to the SdS value, jppn = 1, than the 
Hu-Sawicki mean value is just inside the surface. Mean 
Ippn value for the Hu-Sawicki polytropes converges to a 
value, that lies within the observational bounds, but is 
always different than the ACDM value. 

We found that the mass of a f{R) polytrope is al- 
ways closer to the observed solar mass than the mass of 
the Eddingon polytrope. The difference in the masses 
M H s/Mq ~ 1.002 and M GR /M G ~ 1.048 is systematic 
for all the n p = 3 polytropes. Although the solutions are 
similar with respect to the initial values and obtained 
radii and masses, the solutions never identical. 

The /(i?)-model is stabile with respect to short- 
timescale instabilities [JD] if f"(R) is positive. For the 
present model, this stability condition is satisfied with ex- 



FIG. 3: Approximately one third of the polytropic solutions 
lie within the Cassini limits (j PPN -l)= (2.1 ± 2.3) x 1CT 5 , 
shown in FIGj2] The data concentrates close around the gen- 
eral relativistic value 7ppjv=1, which is a special feature of 
this model. 





ACDM 


f(R)HS 




R/R Q M/Mq 


R/Rq M/Mq -(ppn - 1 


0.5 
0.6 
1.0 
1.2 
2.0 


2.030 1.0484 
1.692 1.0484 
1.015 1.0484 
0.846 1.0484 
0.508 1.0484 


2.000 1.0025 0.00053604 
1.666 1.0025 0.00001434 
1.000 1.0024 0.00040438 
0.833 1.0024 0.00233362 
0.500 1.0024 -0.00005374 



TABLE II: GR values are compared here to modified gravity 
/(-R)HS-model v alue s for different fiducial central densities 9q, 
parametrized in ( |17| |. f(R) polytropes with 60 — 0.5 and 2.0 
are also shown in Fig. |T]). ACDM values were computed for 
reference with the same code. The Lane-Emden solution is 
very similar to 9o = 1.0 ACDM solution in the table. Also, the 
spacetime parameter is not dependent on the initial density 
parameter, bacause of the instabilities in the curvature near 
the center. Note, the space-time parameter for ACDM model 
is (7ppjvT)=— 8 x 10~ 7 in our work. 



tremely small f"(R) everywhere inside and near the poly- 
trope. The mass term within the corresponding Brans- 
Dicke framework, is dependent on a l/f"(R) term. The 
mass of the scalar, m^., is allways big enough for the 
"fifth force" problem to be evaded for these polytropic 
configurations. 



V. CONCLUSIONS AND DISCUSSION 

Although, Hu-Sawicki polytropes are very similar to 
the young sun like stars in GR, only one third of the 
HS polytropes in a typical 29 000 polytrope run have 
7p pn parameter within the experimental Cassini limits 
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Many f(R) models result in 7 PPA r = 1/2 OES], far 
from the GR value jppn = 1- The studied HS-model 
is different in this respect although not fully consistent 
with the observations. For the Hu-Sawicki -model, 'Jppn 
value is nearly one everywhere in our model parameter 
space. However, the f(R)ns space-time around a spheri- 
cally symmetric configuration with a polytropic equation 
of state nearly always converges to a different space-time 
than is allowed by GR. The non-linearity of the system 
poses difficulties, and very high accuracy was not feasible. 
With a sufficient numerical accuracy all the n < 9 models 
were solvable. Furthermore, even as the numerical accu- 
racy was increased JppNjHS value does not converge to 

lPPN,ACDM- 

With only 40% of the jppn values reaching the obser- 
vational limits with the Hu-Sawicki gravity, it is obvious 
that HS-model does not easily reach the Schwarzchild-de 
Sitter solution outside a polytropic stellar configuration. 

This model, therefore, seems to require some fine- 
tuning of the initial scalar curvature in order to repro- 
duce general relativistic SdS space-time outside the sur- 
face of a sun-like HS polytrope. The scalar curvature 
is always unstabile near the center, so no specific value 
can be attached to the HS polytrope, that would exactly 
correspond to the general relativistic polytropic sun. 

To get away from this fine-tuning a physical super se- 
lection rule might exist. We do not suggest any mech- 
anism for this, however, we would like to note that the 
sun is the only star for which the SdS space-time is mea- 



sured to be correct with good accuracy. Could it be, 
that other sun like young stars that can be modeled with 
a polytropic equation of state, might arrive at a slightly 
different vacuum than SdS? In this case HS -model would, 
indeed, describe the space-time around young main se- 
quence stars more accurately. 

The short-timescale instabilities are avoided with ex- 
tremely small and positive f"(R). Also, the scalar field 
mass within the Brans-Dicke scenario becomes large, as 
m a contains a term proportional to 1/ ' j"(R). 

The solutions do not change notably in the mass, ra- 
dius and jppn when the Hu-Sawicki parameters are var- 
ied. This indicates that only the ratio c\/c2, which is 
fixed, is effective for R^> m 2 . However, the "fppN value 
is naturally dependent on the initial central density po, 
that was varied throughout the numerical studies. 

With a different polytropic equation of state, some 
white dwarfs and neutron stars can be modeled. Work is 
under way to test the space-time specifically around f(R) 
low density white dwarfs and neutron star -polytropes. 
Moreover it would be interesting to study numerically 
also other proposed well-behaving action functions. 
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